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ABSTRACT 

This report describes the procedure utilized in the assembly and 
evaluation of a small intermittent blowdown supersonic wind tunnel. 

The tunnel proper was fully constructed and nozzle blocks for 
attainment of a Mach number of 1.5 or 2.37 in the 1/2" by 1.718" 
test section were available. The basic schlieren system was purchased 
complete. An attempt was made to use or convert readily available 
materials for all components. 

It was expected the results of this project would be of a perman- 
ent nature and would be utilized by the Mechanical Engineering Depart- 
ment as a laboratory 1 WeCanietian The installation is suitable for 
the study of compressible fluid flow and the phenonemon of normal 
shocks but is not suitable for the study of the performance of shapes 
or models in the fluid flow. 

The tunnel installation operates satisfactorily with a maximum 
Mach number of 2.29 with approximately 0.28 lbm of air per second. 
Schlieren photographs of the flow field may be taken in Gutter black 


and white or in color. 
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INTRODUCTION 

There are two basic types of wind tunnels available for the study 

of compressible fluid flow (1) the continuous and (2) the intermittent, 
[11] *, With continuous wind tunnels the en of matching com- 
pressor, tunnel volume flow and pressure ratio characteristics present 
major difficulties in the design work. In addition the Mach number span 
of a continuous wind tunnel is usually small and does not extend from 
low subsonic velocities to high supersonic velocities. These factors 
combined with a higher cost of installation and operation of the con- 
tinuous tunnel have resulted in a shift towards the use of the inter- 
mittent tunnels. . 

The intermittent tunnel is further divided into two basic types'*(1) 
the vacuum storage drive and (2) the pressure storage drive, [11] . Gn 
the first system, atmospheric air passes through a drier, wind tunnel 
and diffuses into a vacuum vessel which is continuously evacuated to the 
atmosphere. In the second system, dry air is stored under pressure in 
a storage tank; the air passes cnEeten & throttle valve, wind tunnel. 
and exhausts through a tunnel diffuser to the atmosphere. 

It was the objective of this thesis to install and evaluate the per- 
formance of a small pressure storage type supersonic wind tunnel system. 
The components available at the start of the project consisted of: high 
pressure compressor; storage tanks; two regulator valves; tunnel test 
section; and Schlieren system components; all other pieces for the 
system had to be either specified or designed and specially constructed. 


It was further necessary to carry out assembly and installation of these 


components. Having completed this phase of the objective, it was then 


* Numbers in brackets refer to the Bibliography on page 20. 
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possible to undertake a performance evaluation program. This con- 
sisted of both a qualitative and an quantative approach, The qualita- 
tive evaluation was derived primarily from Schlieren photographs, while 
the quantitative results came primarily from pressure surveys. 

Figure 1 is a schematic diagram of the system showing the basic 
design and necessary components. Discussion of the design and pere- 
tinent details of certain sections of the system is found in Part I. 

Figure 2 is a subassembly drawing of the small tunnel test section 
for which the performance tests were made. Results of these tests are 
presented in photographs, graphical form, and tabular form. Discussion 
of these tests and results appear in Part II. 

The Appendicies include various sample calculations, a discussion 
of Schlieren systems, operational instructions, design of a differen- 
tial pressure transducer, and design of needed heat capacitors for the 


air receivers. 





Part I. Tunnel System Design 
A. Air Supply. 

Air for the system was supplied by a Chicago Pneumatic two stage 
inner cooled air compressor rated at 114 cubic feet of free air per min- 
ute or 0.142 pounds per second. The compressor is driven by a 50 horse- 
power induction electric motor. The compressor discharge pressure can 
be controlled up to 500 psi by the regulator valve (Bailey Type-XSDS-9M). 
Operational instructions for the compressor are given in Appendix A. 

Air from the compressor goes to an after-cooler (Graham Heliflow. 
Type 9AA 145) where the temperature is reduced to about 60°F, and then 
to a cyclone type separator with trap (Wright Austin), where the con- 
densed water is removed. Calculations show the moisure content after 
the air leaves the separator to be about 0.000319 lb. water per lb of 
dry air, Appendix B. Drying of the air is necessary to prevent condensa- 
tion from occuring in the tunnel test section, [6 } ; (11 } . [12] : 

From the separator, the air passes through an oil filter (C. M. 
Kemp, Type FPF-1610), using activated alumina (Vapoilsorb) as a desi- 
ccant. This filter was modified by putting a 60 mesh stainless steel 
strainer on the inlet side along with a fiberpacked strainer, This modi- 
fication was made to prevent blowback of desiccant particles to the 
compressor by high pressure air when the compressor was unloaded. 

The air then goes through the regulator valve, which maintains 
a constant compressor discharge pressure, to the two air receivers, 

‘The receivers, connected in parallel, provide a total storage volume of 
234 cubic feet. Figure 3 is a photograph showing these various compon- 


ents in relation to the acoustic enclosure which houses the compressor. 
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When operating the wind tunnel, the mass flow of air from the 
receivers to the tunnel section is controlled by a regulator valve 
(Bailey, Type XSDS-9M). The regulator valve provides for stagnation 
pressure control up to 200 psi. Figure 4 shows the position and loca- 
tion of this valve. 

B. Tunnel Section. 

The tunnel section consists of six basic parts: stilling tank; 
entrance nozzle; test section; diffuser; flow meter; and back pressure 
control. 

The stilling tank was made from a 30 inch piece of nominal 6 inch 
pipe flanged at both ends, The flanges were faced and grooved for use 
with O-ring seals, The upstream end of the tank was connected by 
means of a piece of high pressure rubber hose to the regulator valve 
discharge. Figure 5 is a photograph of these pieces, In the upstream 
end of the stilling tank, 3 stainless steel 60 mesh screens were used 
to diffuse any high velocity profile. 

The entrance nozzle was constructed of a polyester resin and 
chopped fiber glass molded over a wooden core. A pair of cubic equa- 
tions were used to define the nozzle contour so that flow separation 
would not occur. The flanges were molded integral with the nozzle and 
were faced and drilled after casting. The upstream flange has a 1/2 
inch lip for a matching connection with the stilling tank. 

The test section was designed and constructed by the Ames 
Aeronautical Laboratory, Moffett Field, California. Features of this 
test section consist of a variable area supersonic diffuser, ‘aban 
changeable nozzle blocks, and interchangeable side walls for either 


visual or pressure surveys. 
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The present sets of nozzle blocks are for Mach numbers 1.50 and 
2.37. The sidewall for pressure surveys is provided with 32 pressure 
taps. Figure 2 is a subassembly drawing of the section and Figure 6 
is a photograph showing the section with side walls removed, Tables {I 
and II give the nozzle coordinates, area ratios, and design Mach numbers 
and pressure ratios based on one-dimensional calculations for the two 
sets of nozzle blocks. 

The diffuser, like the entrance nozzle, was also molded over a 
wooden core with a polyester resin and fiber glass. To provide added 
strength, 1/8 by 3/8 inch steel bars were molded into the side and 
flanges. The design was strafght wall with a maximum divergance angle 
of 7° to prevent flow separation; and diffusion of the air stream was 
carried to a point where the flow could be treated as incompressible for 

_metering purposes. 

The four pieces just described were integrally mounted to a movable 
stand and are shown in Figure 7. 

The flow meter was designed according to A.S.M.E. standards for a 
sharp edge orifice with D and (1/2)D pressure taps, [4] 5 [5] . An en- 
trance length of 20 pipe diameters and an exit length of 10 pipe diameters 
were provided. Nominal 3 inch pipe was used which was bored to a 3.120 
inch diameter for 16 inches either side of the orifice. O-ring seals 
were used throughout including the flanged joint with the diffuser. 
Orifice plates were made of 1/8 inch stainless steel. Appendix C presents 
the details for computing mass flow rates, 

The back pressure control consisted of a manually operated 3 inch 


globe valve. This control was incorporated so that changes could be made 
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in the test section air density. Since this throttling was somewhat 
noisy, a large truck muffler was added to the valve exit. 

These pieces were also mounted on a movable stand and are shown 
in Figure 8. 

The completely assembled tunnel section is shown in the photo- 
graph, Figure 9. 
C. Instrumentation. 

Instrumentation for the tunnel section consisted of devices for 
three types of measurement; pressure, temperature, and velocity. 

Manometers were used for all quantitative pressure measurements 
while bourdon tube gages were used for qualitative check purposes. 
Two 120 inch mercury differential manometers were utilized with the 
32 tap pressure survey sidewall of the test section. Two 16 connector 
manifolds with quick opening valves were made to accommodate the 32 taps 
to the two manometers. Figure 10 shows the arrangement of the manifolds 
with the manometers. It was originally planned to produce strain gage 
type differential pressure transducers for these measurements to de- 
crease the reading time required; fabrication difficulties, however, 
prevented the completion of this project, Appendix D. 

A 120 inch mercury differential manometer was also used to monitor 
the test section throat pressure continuously during pressure survey runs. 

Stagnation pressure, from the stilling tank tap, was likewise con- 
tinuously monitored with a 120 inch mercury manometer besides being fed 
to the differential manometers. 

Static pressure for the flow meter could be read with either an 
oil (sp. gr. = 0.834) or mercury filled 60 inch manometer, Differen- 


tial pressure for the orifice was determined with a water manometer, 





Ellison inclined manometers were used in conjunction with the pitot 
tube, 

The only quantitative temperature measured was the stagnation 
temperature in the stilling tank. This was accomplished with a 
stainless steel radiation shielded 30 gage copper-constantan thermocouple 
referenced to ice. The output emf of the thermocouple was read with a 
continuous balance precision potentiometer (Minneapolis-Honeywell, Model 
SY156X-VH). Various other temperatures of the tunnel system were meas- 
ured with gas bulb indicating thermometers. 

Velocity measurements in the stilling tank were made with a 1/8 
inch diameter Prandtl-Pitot probe. These meaSurements were made sole- 
ly for the purpose of determining the uniformity of the velocity pro- 
file leaving the straightening screens in the stilling tank and were 
only of a qualitative nature. The Pitot probe was removed from the 
stilling tank after finding that the screens adequately diffused the 
high velocity entering air to a uniform low velocity profile. 

Velocity measurements and visual flow studies in the test sec- 
tion were accomplished with a coincidence type Schlieren systen manu- 
factured by Tinsley Laboratories Inc. Aschematic diagram is shown in 
Figure 11. The components were mounted to a rigid moveable stand with 
the wheels set in tracks parallel to the test section, Figure 8. The 
viewing section of the Schlieren system was modified by installation 
of an Ilex Synchro shutter and a 4 by 5 inch ground glass view plate 
and film pack holder, Figure 12. This modification allowed photographs 
of the flow behavior to be made which, when enlarged in printing, pro- 


vided for qualitative and some quantitative measurements. Appendix E 





presents a general theory of the Schlieren method of flow visualization 


as well as operating instructions for this particular system. 





Part II. Performance Evaluation 


A. Qualitative Results. 


~ 


The qualitative evaluation of the tunnel system performance was 
derived primarily from Schlieren photographs. Figure 13 depicts the 
tunnel section operation under shock free conditions. The knife edge 
intercepted the light beam in both a vertical and horizontal position 
to show the change of density patterns in the horizontal and vertical 
directions, Figure 13c, with the knife edge horizontal, clearly shows 
the boundary layer. 

Criss-crossing the tunnel outline are Mach wavelets produced by 
wall friction and other irregularities of the tunnel surfaces. These 
wavelets are related (2] to the Mach number of the air flow in the 
tunnel by the relationship, 

sin & = be 

M. 
where X is the acute angle between the wavelets and a horizontal line. 
Figure 18a is an enlargement of a section of Figure 13a and Figure 18b 
is an enlargement of Figure 13b. The angle G in Figure 18a is 27.50° 
with a resulting Mach number of 2.17. This wavelet originated at the 
junction of the nozzle and the test section. Utilizing the pressure 
ratios, the Mach number at the corresponding point was determined to 
be 2.125. The angle in Figure 18b is 29,25° with a resulting Mach 
number of 2.04. The Mach number calculated for this position was 1.828. 
This wavelet originated at approximately pressure tap 30. 

Figure 14 depicts a normal shock occuring in the diffuser section 
which was being held at a constant height of 1.718 inches. Inner action 


between the shock and the boundary layer is clearly indicated. Up- 


Stream of the shock the boundary layer is stable. At station 12 there 
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is a separation of the boundary layer due to the interaction of the 
shock and the boundary layer. Less severe interaction occurs between 
the boundary layer and the repeated shocks downstream until complete 
subsonic flow is established downstream of station L5. 

Comparing Figure 14a operating under a stagnation pressure of 
67 inches of mercury and a back pressure of 30.5 inches of mercury to 
Figure 14c for tunnel operation with a stagnation pressure of 81.5 inches 
of mercury and a back pressure of 38.3 inches of mercury, it may be seen 
that the resultant shock has moved upstream almost three inches. The 
pressure ratio at station 11 in Figure 14a where no shock exists is 
approximately 0.1088 as taken from Figure 23 or Table V. At a point 
forward of the shock in Figure 14c a pressure ratio of 0.103 exists 
as taken from Table VII.The increased density of the air and the high- 
er Reynolds number must account for the shock moving upstream. 

Reduction of the stagnation pressure also moves the shock up- 
stream as shown in the sequence of photographs in Figure 15. A strong 
shock exists at the nozzle exit of Figure 15a with a stagnation pressure 
of 65.35 inches of mercury. Reduction of the stagnation pressure to 
59.96 inches of mercury allowed the shock to move upstream one and one 
half inches towards the nozzle throat as shown in Figure 15c. Reduction 
of the shock intensity may be noted also as the shock moves upStream by 
comparing Figures 15a and 15c. 

~The minimum starting area of the variable diffuser was found to 
range between 0.695 and 0.700 square inches as compared to the isentropic 
minimum area of 0.665 Square inches, The area actually required corres- 
ponds to a diffuser heights of 1.39 to 1.40 inches. Figure l6a is a 


photograph of the diffuser nozzle set at a height of 1.39 inches. The 
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shock could not be blown through the diffuser with a stagnation pressure 
of 82.71 inches of mercury. Increasing the diffuser height to 1.40 
inches, (area change of 0.005 square inches) permitted the shock to be 
blown through the diffuser nozzle and the stagnation pressure reduced 


as shown in Figure 16b. The stagnation pressure could have been reduced 


to move the shock farther upstream to the throat position. Figure 17a 
and 17b depict tunnel operation with the diffuser height set slightly 
greater than the required minimum and show the shock standing in the 
diffuser throat and the stagnation pressure reduced to a minimum while 
maintaining the shock location. Several minutes after Figure 17a was 
taken the stagnation pressure was reduced from 60.00 inches of mercury to 
59.50 inches of mercury and the shock moved through the diffuser nozzle 
and relocated itself in the diverging portion of the tunnel nozzle. 

The sharpness of the Schlieren photographs was poor due to the use 
of non-collimated light and the fact the light beam passed through the 
test region twice, Appendix E. Figure l6c gives an indication of the 
seriousness of this problem. The knife edge was not inserted into the 
light beam. The black parallel vertical lines are images produced by 
placing fine copper wires on the outside of the Schlieren window nearest 
to the spherical mirror; the diagonal parallel lines are images produced 
by placing the same size wire on the outside of the Schlieren window 
nearest to the light source. The greater the distance between the 
images of any wire, the farther the light beam had to travel before 
reaching the spherical mirror. The non-parallel light beam intercepted 
the diagonal wire first and hence the image of this wire was deflected 


more due to the longer distance to the spherical mirror than it was for 
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any of the vertical wires which were closer to the spherical mirror. As 
evidenced by the two images formed for each wire, the Schlieren picture 
of the of the air flow would tend to be fuzzy due to superimposing the 
two images formed. The horizontal black line is a wire placed at the 
tunnel centerline. 

Also contributing to poor quality of the Schlieren pictures was 
the incomplete removal of oil vapor from the air. During tunnel opera- 

\ 
tion, oil would be deposited in a very thin film on the glass wall and 
cause a refraction of the light. In addition, the film deposited was 
not of a uniform nature. This condition may be seen in Figure l7c. 
The dark convex lines at et at Mbie 4 and 5 are the result of an oil film 
deposit. These same lines could still be seen after securing the air 
supply and leaving the Schlieren unit in operation. Due to this oil 
film formation, it was necessary to remove the Schlieren window at 
frequent intervals and wipe off the small amount of oil that had been 
deposited. Water vapor in the atmosphere condensing on the outside 
surface of the Schlieren window was also a factor in obtaining clear 
Schlieren photographs. 

As the vertical wires were spaced at one inch intervals, the dis- 
tance between wires in the photograph may be measured and comparison of 
this measured value to the known one inch dimension permits determining 
a scale factor for the photograph. This method was utilized to deter- 
mine the diffuser heights for heights other than 1.718 inches. A 


lucite gage block, 1.718 inches high, was utilized when runs of this 


diffuser height were made. 
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The apparent difference in intensity of the shock and 
the visibility of the boundary layer caused by the knife 
edge orientation demonstrates why a horizontal knife is 
generally used when studying boundary layers and a vertical 


knife is used when studying shock waves. 


B. Quantitative Results. 

The quantitative evaluation was derived primarily from 
pressure surveys. The Mach numbers actually obtained in 
=a tunmel were generally less than the predicted adiabatic 
Mach numbers due to friction effects not accounted for, ob- 
lique shocks and other departures from a one dimensional 
flow analysis. Figure 20 is a plot of the adiabatic Mach 
numbers obtained from the pressure and area ratios for a 
constant diffuser height of 1.718 inches vs. distance down- 
stream of the tunnel nozzle throat. The solid curve is the 
isentropic Mach number based on the change of area ratios 
only as given in Table I, while the dashed curve is the 
predicted adiabatic curve as determined by the use of iso- 
Cclines in Appendix F. The circles are the adiabatic Mach 
numbers obtained from the data taken on the tunnel and are 
tabulated in Table V, based on the procedure outlined in 
Appendix F. 

Up to a Mach number of 1.936 at pressure tap number 
4, the actual Mach numbers followed the adiabatic curve 
closely, exceeding both the isentropic and adiabatic 
curves between pressure taps numbers 4 and 7. At this 


point, the Mach numbers obtained in the tunnel began to 
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fall below the isentropic curve but still above the adia- 
batic curve to pressure tap number 9. The maximum Mach 
number of 2.275 was obtained at pressure tap number 7. The 
decrease in Mach number at taps number 10 and 24 is un- 
explained. This decrease in Mach number at these pressure 
taps appeared consistently in all runs made during the 
evaluation of the tunnel. 

An attempt was made to maintain a constant Mach number 
in the test section by increasing the height of the varia- 
ble area diffuser. The results of this run are shown 
graphically in Figure 21 and in tabular form in Table VII. 
Because of design limitations, the diffuser height could 
not exceed 1.883 inches. The calculated maximum height of 
the diffuser to maintain a constant Mach number in the test 
section including the effects of friction was 2.145 inches, 
Appendix F, based on a friction factor of Afie=sOeOle Table 
IV presents the diffuser coordinates vs. diffuser length in 
tabular form while Figure 19b presents the same data in 
graphical form. Examination of Figure 21 indicates a 
definite drop in Mach number between the nozzle exit and 
tap 24, in much the same manner as the drop for constant 
diffuser area, Figure 20, between the same pressure tap 
locations. As the diffuser did not increase in height 
appreciably in the first 1.5 inches of diffuser length, the 
effect of friction was to continually reduce the Mach num- 
ber. Recovery of the Mach number was accomplished shortly 
after the diffuser height began to increase. The recovered 


Mach number is somewhat higher than that for constant 
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diffuser area, and appears to be maintaining a more con- 
stant value, thus indicating the effect of increasing the 
area of the diffuser. Had it been possible to increase the 
diffuser height to the calculated height of 2.145 inches, 
the resulting Mach number would have been increased much 
more. Until the behavior in Mach number without area 
Change can be explained, the effect of changing the area of 
the diffuser can not be adequately predicted. 

Figure 24 is a plot of adiabatic Mach numbers vs. 
distance from the tunnel nozzle throat for two operating 
conditions: (1) the diffuser height set at constant 1.718 
inches and a back pressure of 38.98 inches of mercury, and 
(2) the diffuser set at the minimum starting height of 1.40 
inches with a back pressure of 30.54 inches of mercury. 

The minimum diffuser height was based on the actual minimum 
height required to operate the tunnel supersonically in the 
test section confirmed by the calculated height in Appendix 
F. The calculated minimum height for isentropic flow was 
found to be 1.33 inches; while that for adiabatic flow was 
1.40 inches, Appendix F. Table III gives the diffuser co- 
ordinates in tabular form and Figure 19b presents the same 
information graphically. It should be noted there was 

very jie enka change of area of the diffuser near the nozzle 
outlet when adjustment of the diffuser height was made 
downstream. 

The data followed the predicted curves in much the 


same general pattern as when the tunnel was operating under 
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a constant diffuser height shock free condition. The 
characteristic drop in Mach numbers may be noted at pres- 
Sure taps 10 and 24. A predicted adiabatic Mach number at 
the diffuser throat was calculated with the diffuser height 
set at 1.40 inches, Appendix F. This predicted value of 

M = 1.64 does not correlate too well with the value of 

M = 1.34 obtained in the tunnel as shown in Figure 24. 

While operating supersonically, the tunnel had a mass 
flow rate ranging from 0.25 lbm/sec. to 0.36 lbm/sec., as 
determined by the use of an orifice meter. The correspond- 
ing stagnation pressure varied from 62 to 88 inches of 
mercury. The mass rate of flow for the orifice meter may be 
determined by the procedure outlined in Appendix C, or the 
imal nass rate of flow for the tunnel nozzte"throatyas 
determined by Fliegner's formula, may be corrected by the 
wee of the»empirical Coefficient of Dirscharge presented 
graphically in Figure 25. 

The average friction factor 4f, in the diffuser sec- 
tion was found to be 0.01865 when the tunnel was operating 
shock free supersonically, Appendix F. The diffuser height 
was maintained at a constant 1.718 inches and the data 
from Table V utilized. All predictions made involving the 
use of friction were based on the assumption the mean 
friction coefficient 4f was 0.01 and was constant at this 
value. 

Figures 22 and 23,sand Table» ixXepresent the resul@e 


of a pressure traverse of the tunnel operating with a 
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normal shock in the tunnel nozzle. Figure 22 is a plot of 
the Mach number vs. distance from the nozzle throat. The 
solid line represents the isentropic case assuming a single 
normal shock occurring at M = 2.05 and the dashed line the 
adiabatic case assuming a normal shock at M = 2.00. Imme- 
diately downstream of the single shock the isentropic Mach 
number will be decreased at 0.569 and continue to decrease 
in the nozzle, due to the increasing area, until it reaches 
a value of 0.388 at the end of the nozzle; the sen taepic 
Mach number then remains constant in the diffuser section. 
For the adiabatic case, the Mach number immediately down- 
stream of the shock will be 0.577, based on the shock occur- 
ring at M2» 2.00. The Mach number continues to decrease for 
the subsonic flow but at a slower rate. Between pressure 
taps 9 and 15, the area change is slight and the Mach number 
remains constant at M = 0.429. In the constant area dif- 
fuser, the adiabatic Mach number is increased due to the 
frictional effects. As shown in the Schlieren photographs 
by the repeated shock pattern, Figures 13 to 18, the flow 
conditions through a normal shock are definitely not a 
Simple one dimensional process. A shock occurs and the flow 
immediately goes subsonic. Due to inner action with the 
boundary layer, the flow accelerates to a lower supersonic 
value and the process iS repeated until the flow remains 
Subsonic. The overall effect of Mach number reduction 
caused by the multiple shocks is shown in Figure 22. Due 


to the complex inner action of the normal shocks and the 
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Oblique shocks with the boundary layer, and the change of 
stagnation pressure for each shock occurring, definite 
determination of the Mach number through a shocking system 
is difficult. The method for calculating the Mach numbers 
of Figure 22 is presented in Appendix F. 

At a Mach number exceeding 1.00 the effect of friction 
ESeEO seduce, the Mach number to,a limiting value of 1.00. 
At Mach numbers lower than unity, the effect of friction is 
reversed; namely, the subsonic flow is accelerated to a 
limiting Mach number of 1.00. This increase of Mach number 
Gue@eto the effects of friction on subsonic flow may be noted 
by the gradual increase of Mach number of the subsonic flow 
in Figure 22 or in Table IX. 

Figure 23 is a plot of the tunnel pressure ratios 
based on the stilling tank stagnation pressure vs. distance 
from the nozzle throat for the cases of shock free and 
shock in the nozzle operations. The pressure ratio pattern 
for shock free operation follows the isentropic curve 
closely with deviations commencing about four inches from 
the nozzle throat. With shock in the nozzle operation, 
failure of the pressure ratio downstream of the shock to 
recover to the isentropic value is a loss of stagnation 


pressure caused by the repeated shock pattern. 
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Conclusions and Reconnferadtadhs 
A. Conclusions. 
1. The qualitative and quantitative evaluation of the tun- 
nel system as installed, show it to be a suitable laboratory 
device for the study of compressible fluid flow. 
2. The minimum starting height of the variable area dif- 
fuser was found to be 1.40 inches. 
3. The average friction factor, 4f, was found to be 0.02044, 
but is believed high in value. 
4. Schlieren photographs of the flow field may be taken in 
black and white or in color. 
5. A vertical knife edge orientation in the Schlieren sys- 
tem is best for studying shock phenomenon. 
6. A horizontal knife edge orientation in the Schlieren 
system is best for studying the boundary layer. 
7. #The sharpness of the Schlieren photographs is poor due 
to the use of non-collimated light and the light beam pass- 
ing through the tunnel test section twice. 
8. The phenomenon of a standing shock is definitely not 
one dimensional. 
B. Recommendations. 

It is recommended the following modifications for im- 
proving the performance and usefulness of the tunnel system 
be undertaken: 

1. An additional relay be installed on the Bailey stagna- 
tion pressure controller to maintain a more uniform stag- 


nation pressure. 
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2. Additional heat capacitors be installed in the receiv- 
ers to maintain a more constant stagnation temperature. 

3. A total pressure probe be installed in the fixed geom- 
etry diffuser to obtain the back pressure stagnation pressure. 
4. For more quantitative visual results, the Schlieren 
system be modified to a single pass parallel light beam 
system. 

5. A more rapid means of obtaining the differential pres- 
sure measurements be installed. 

6. The camera system be modified to a reflex type system. 
7. An additional oil filter be installed. 

8. Additional pressure taps be installed downstream of the 
variable area diffuser in the existing pressure survey tun- 
nel sidewall. 

9. The geometry of the variable area diffuser be modified 


to permit constant Mach number runs. 
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Figure 5. Stilling Tank. 





Figure 10. Manometer Manifolds 
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Figure 12. Schlieren System Light Source, Knife Edge, 
and Camera-Viewing Components. 
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A. Vertical Knife. Po = 71.46 in Hg. @., Ph = 30.92 in, Hg. @., Ty = 61.1° F. 
Nozzle Running Shock Frea with Air Flow Right to Left. Nota Strong Mach 
Wavalat Originating at Junctinn of Nozzle and Tast Sactinn. Saa Figura 18. 





B. Vertical Knife. P, © 72.31 in Hg. a., Py © 30.92 in. Hg. &., To © $9,23° F. 
Tunnel Running Shock Free With Air Plow Right to Left. 
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C. Horizontal Knife. P, © 72.71 in. Ag. 8., Py © 30.92 in. Hg. @., Ty © 65.5° F. 
Tunnel Running Shock Fraa With Air Flow Right tn Left. Mota Thicknasa nf the 
Boundary Layers. 


Figure 13, Tunnel Running Shnck Pras 
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A. Vertical Knife. P, = 67.35 in. Hg. @., Pp = 30.45 in. Hg. a., 7, = 6653° F. 
Diffuser Height 1.718 inchee, Nnrmal Shock et Stetinn 12. 
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B. Horizontel Raffe. P, «© 67.41 in. Ag. @., Pp = 30.86 in. Ag. @., T, = 65.9° F. 
Diffueer Height 1.918 inchee, Normal Shock at Station 12. 





C. Vertical Knife, ra 81.52 in. Hg. @., Ph = 38.32 in. Hg. a., T, = 57.0° P. 
Diffueer Height 1.718 inchee. Due to Higher Deneity and Reynolds Number 
Shock hae Movad Upetraam ae Compared tn Figura 14A, 


Figure 14, Effect of Increased Back Preseure nn Shock Locetinon 
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A. Verticel Knife, P, © 65,35 in. Hg. a., Pp © 30.65 in. Hg. @., Ty = 59.6° F. 
Boundery Layer Vieible on Nnzzle Walle end Incraesing in Thicknese Downstraen 
of Shnck,. 





B. Horizontel Knife. Po © 64.41 in, Hg. a., Py © 30,83 in, Rg. a., To © 63.2° F. 
Note Pull off of Boundery Layer et Shock end Asyometrical Shading of Cnrras- 
ponding Perts on Oppneite Sides of Tunnel Centarline. Shock ie Moving Tnvarde 
Nozzla Throet es P, is Decreesad from Figure 15 A. 





C. Vertical Knife. P, = 59.96 fo. Hg. a., Py, - 30.73 in. Hg. a., T, = 61.8° F. 


Weaker Shock Moving Upstreem es P, is Reduced en Additionel Amount. 


Figura 15. Effect nf Reducing Stagnatioo Pressure for Shnck io Nngzle 
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A. Verticel Knife. P, = 82.71 in. Hg.e., Ty = 57.6° F, Diffueer Height © 1.39 in. 
Shock Could not be Blown Through Diffuser Nozzle. 





Verticel Knife. Po = 60.26 in, Hg. @., Ty = 58.6° F. Diffoser Height = 1.40 fo. 
Aree hae been Increesed 0.005 Square Ioches Compered to Figure 16 A, eod Shock 


hse been Blown Through Diffuser Nozzle. Pressore Coold be Reduced more to move 
Shock towards the Throat. 





C. Tunnel oot Running. Knife Edge not Ioserted. Effect of Noo-Perellel Two Pess 
Light ie Evidaot. Double Verticel Lices ere Images of Wires Pleced oo 
Sphericel Mirror Side of,Schlieren Wiodow, Diegonel Liocesa ere Images of Wiree 
Pleced on Light Source Side of Schlieren Window. Light Pessing These Wires 


hed to Trevel Ferther before -Striking Sphericel Mirror end Hence More Deflection 
of the Imace. 7 


Figure 16. Effect of Aree Change on Shock Location in Diffuser Nozzle, 
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A. Vertical Enife. P, = 60.00 in. Ag. @., To = 55.0° F. Diffuser 
Reight © 1.42 inches. Shock Pulled Back Through Diffuser Nozzle 
wheo P, was Decreased to 59.30 in. Eg. e. 





B. Borisontel Mnife. P, © 60.61 fo. Ag. a., Ty © 63.1° F. Diffuser 


Height © 1.47 fechas, Rote Strong Pull Off of Boundary Layer et 
Station 12.25. 





C. Vertical Knife. P, © 72,31 in, Hg. Ge, T> oe $8.7° F. 
Lines at Statiooa e 


Dark Horizonte! Lines Sterting et Stetion 3 ere elso off Stresks on 
Scbhlisres Welle. 


Dark Concave 
end 5 ere the Result of O11 Film Deposit oo Walle, 


Figure 17, ‘Twnmel Renming With Redeced Diffuser Area. 
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A. Verticel Knife. P, = 71.46 io. Hg. e. Py = 30.92 fo. Hg. @. T= 61.1° FP, 


Air Plow ie Right to Left. Section ie eo Lolergement of Figure 13 A, 
Nozzle is Running Shock Free. Mach Number determioed from Wevelets by 
Ralaticosbip Sina = 45 « At Corraepooding Potot Adiebetic Mach 
Number ie 2.125. 





Varticel Koife. P, = 72.31 10. Hg. @., Py = 30.92 fo. Hg. @., Ty = 59.23° P. 
Air Plow ie Right to Left. Section te an Zolergeweot of Figure 13B. At 
Correepooding Point Adiabatic Mach Mumber is 1.828. 


Figure 18, Determinatioo of Mach Number from Mach Wevelete,. 
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Vertical Knife. P, = 71.46 to. Hg. a. Py = 30.92 to. Hg. @. T= 61,1° F, 
Air Flow fe Right to Left. Saction is eo Zolergewent of Figure 13 A, 
Nozela ie Running Shock Free, Mach Number detarmined from Wevalete by 
Reletfonship siq x = 4% « At Corresponding Point Adfebatic Mach 
Number fe 2,125. 





Varticel Knife. Po = 72.31 in. Hg. e., PL = 30.92 in. Hg. Mey To = 59.23° | 
Air Flow fe Right to Left. Sectfioo is an Enlergament of Figura 13B. At 
Corrasponding Point Adiabatic Mach Number fe 1,828. 


Figure 18, Determinetion of Mach Numbar from Mach Wevelets. 
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APPENDIX A 


TUNNEL SYSTEM OPERATION PROCEDURE 


I. For charging the storage tanks the following procedure should be 


followed: 


_ 


10. 


Lae 


ne 


Open the 3/4 inch valve connecting the trap with the separa- 
tor. 
Close the 3/4 inch drain valve located below the tee after 
the trap. 
Open the 1/4 inch balance line valve between the separator 
and the trap. . 
Open the cooting water valve to the after-cooler heat ex- 
changer. 
Close both two inch inlet valves located at the entrance to 
the low pressure storage tanks, 
Open the valve on the low pressure after-cooler. 
Items 5 and 6 are safety precautions to prevent high 
pressure air from entering the low pressure tanks. 

Close the drain valve on each high pressure tank. 
Close the two inch discharge valve on each high pressure 
tank outlet. 
Open the inlet valve to each high pressure tank. 

Place the unloader on the compressor in the START position 
(vertical with the index head down. ) 

Place cover on the compressor oil vent. 

Open the air supply valve to the Bailey controller (Pressure 
should read 28 psi). 
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14. 


at. 


16. 


MT 


18. 


19. 


Zz. 


Zt. 


Adjust compressor discharge pressure for a pressure slight=- 
ly higher than any pressure indicated in the receivers. This 
is done by turning the adjustment knob on the Bailey Con- 
troller counterclockwise until the red hand indicates the de- 
sired pressure, 

Connect the cooling air circulating centrifugal blower in 
the compressor shack, 

Throw the master power switch located in the main bulkhead 
switch panel for the High Pressure Air Compressor to the ON 
position. 

In the compressor shack, throw the overload trip arm to the 
ON position and push the START button. 

After the compressor has started running, place the unloader 
valve in the RUN position (handle horizontal with the index 
head pointing towards the compressor. ) 

Adjust the compressor discharge pressure regulator valve to 
give the desired discharge pressure, (Maximum 500 psi) 

Maintain a check on the storage tank pressure as an automatic 
unloader for the compressor is not provided. 

Should an excess amount of air be coming out of the drain line 
from the trap, it is probable that additional water is required 
in the trap. Water may be added when the compressor is secured 
through the union connection of the balance line until it is 
observed coming out of the trap drain. 


To secure the compressor place the unloader valve in the START 
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position, push the STOP button and remove the oil vent cover. 
An emergency STOP button is also located with the Bailey con- 
troller, 

22. Close the inlet valves to the storage tanks. 

23. Close the balance line valve. 

24, Close the valve connecting the separator and trap. 

25. Close the air supply valve to the Bailey Controller. 

26. Blow down excess air in the lines from the compressor to the 
Storage tank inlet valves by opening the trap drain valve. 

27. Secure the cooling water. 

28. The drain on each high pressure tank and the inner cooler 
should be blown at periodic intervals to remove any accumula- 


tion of oil or water. 
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eee 
ae. 
24, 
Zo. 


26) 


ZB, 


28, 


position, push the STOP button and remove the oil vent cover. 


An emergency STOP button is also located with the Bailey con- 


troller. 


Close the 


Close the 


Close the 


Close the 


Blow down 


inlet valves to the storage tanks. 

balance line valve. 

valve connecting the separator and trap. 
air supply valve to the Bailey Controller. 


excess air in the lines from the compressor to the 


Storage tank inlet valves by opening the trap drain valve. 


Secure the cooling water. 


a 
= 


The drain on each high pressure tank and the inner cooler 


should be blown at periodic intervals to remove any accumula- 


tion of oil or water, 
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Il. 


Blow down procedure 


rt. 


a 


3. 


Open the outlet valves on each high pressure tank. 

Open the air supply valve to the stagnation pressure Bailey 
Controller (Supply pressure should read 28 psi). 

Open the quick opening valve at the outlet of the regulator 
valve. 

Slowly turn the adjusting knob of the Bailey controller 
counter clockwise to allow a gradual buildup of pressure in 
the stilling tank to the desired stagnation pressure, A load- 
ing pressure of at least 5 psi is required before the re- 
gulator valve will open. 

Figure 1-A gives the approximate differential pressure requir- 
ed to obtain the maximum Mach number in the tunnel test sec- 
tion for a given stagnation pressure measured above atmospheric 
pressure in inches of mercury. 

The time of each run _— be prolonged as well as reducing the 
rate of change of stagnation temperature by allowing the com- 
pressor to continue charging the storage tanks while opérat- 


ing the tunnel. 
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APPENDIX B 
AIR MOISTURE CONTENT CALCULATIONS. 


Assume that the perfect gas law relationship applies 


PV = mRT 
where P = gas pressure Ib | Ff 
V = gas volume f= 
m = lb mass of gas 
R =# universal gas constant of gas, Ft - 10g 
ee 
Oo . lO m R 
© 2==qas temperature, Rankine. 


The ratio of mass water vapor in the air to the mass of air 


> 
= 


is then given by 


fe 
My _ fy 7 =- fe 
Me AY Ry, Fa 
fer 


where the subscripts a refer to air and v refers to the 





water vapor. Utilizing the value of R, = 53.34 and 
R, = 3 aa 
| ecules eu Glee 
my = O@rleo 2 = "COL Dp ‘ 
Ng = yqeyr Sed 


Upon exit from the heat exchanger, the air is saturated at 
500 psia and 60°F. From steam tables at 60°F 


Psat = QO.2563 psia 


and 





aa 


4500-02565 | (/b.dryar 





ey = 06225 
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APPENDIX C 
FLOW RATE CALCULATION PROCEDURE. 

The mass rate of flow for the tunnel, as measured at 
the orifice plate, was determined by the procedure out- 
lined by the American Society of Mechanical Engineers [5]. 
All equations, figures or table numbers correspond to num- 
bers in that publication. Due to the removal of the water 
vapor by the air drying system, it is assumed the air is 
ary for metering purposes. In addition, due to the low 


velocity of the air, the air is considered incompressible. 


> 


For a representative run: 


cl = /2.08. in. Gigs 


le 5 oe 25m 


Pob = 30.84 in. Hg. 


D = Pipe diameter = 3.120 inches. 
q = 2.496" = orifice diameter 
_ orifice diameter = d = 0.800 
p ~ pipe diameter D 
AP = differential pressure in inches of HO = 4.55 


Equation 2, page 57, gives the mass rate of flow lbm/hr. 


i 359CF4°F_ Y hw 


49 





Where C = coefficient of discharge 
F — £eleaiuty, Cm ppproach “faccor 
Fe = ternal] gexpdnsion factor 
reers«Cs =~ =s«C«s TT EPerential pressure in THeHes oF oo 
ai )6=0«6 specific weight of ain at the inlet side of the 


primary element, lb/ft 


From Figure 38, page 67, for a stainless steel orifice plate 
at 65°F 
Fa = 1.00 
From Table 9, page 80, at 65°F 
O ade = 0.075650 lbAete at 29.92 in. Hg. 
Correcting for increase in pressure 


Yair = 0.07565 x 32-85 #os0762 leyee 


2 2 


qd = (2.496 ne = Oo. 2350 nil. 


Aesume avalue of K = CF = O.80 and calculate an aGpprox- 


imate m 


3 
il 


359 coe OO x 6.230 V 4755 x O.0762 


1053 lbm/hr. 


And g/ from Figure 35, page 62, for 65°F air = .0000122 


Calculate Reynolds number ne from Equation 8. 


— 0.004244 m _ 0.004244 x 1053 = 117,600 
Bo 3.120 x 0.0000122 


From Table 5, page 26, or from Figure 1-C of this Appendix, 


for Pp 


is K 


0.800, and R = 117,600 the flow coefficient 
D 


0.8062 


Corrected calculation of flow rate, therefore is 


rm 359-x..0.8062 x 6.230 -V 4.55 x 0.0762 


1060 lbm/hr. = 0.2945 lbm/sec. 
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The ideal mass flow rate at the nozzle throat is determined 
by Fliegner's formula [2]. 

° A”P 

“mee 0.582 2 


Where m= ideal flow rate lbm/sec. 

A* = tunnel throat area in Square inches with 
m = 1.00 

Por = stagnation pressure psia 

Be = stagnation temperature in °Rankine 

m: On.532 x 0.367 x 72,508 x (Oem 


‘ 925 


0.301 lbm/sec. 


The nozzle discharge coefficient C, is defined as the 


ratio of the actual nozzle flow to the flow calculated 
from the isentropic laws for the initial and final pres- 
sures of the actual nozzle. Based on this definition the 


coefficient of discharge is 


oe 
c = o8etuet & O.2945 e& 0.962 
P m 0.301 

ideal - 


Figure 1-C presents the Flow Coefficient K plotted 
vs. the Reynolds Number for a flat plate orifice with 
pressure taps at D and 1/2 D for a three inch pipe. This 
figure presents the data of Table 5 (5) graphically for 
values of d/D equal to 0.580 and 0.800. Figure 25 presents 
graphically the nozzle discharge coefficient a plotted vs. 
the stilling tank stagnation pressure. Values of S65 were 
determined by the use of Fliegner's formula and by the 


al 





procedure outlined by the American Society of Mechanical 


Engineers (5) as used in this Appendix. 
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APPENDIX D 
DESIGN OF DIFFERENTIAL PRESSURE TRANSDUCER. 

Differential pressure transducers were originally to have measured 
the differential pressure between the tunnel pressure and the stagna- 
tion pressure. This method of measurement was discarded due to the 
difficulties of fabricating similar output transducers, 

A successful transducer, Figures 2-D and 3-D, was made based on the 
principle of the deflection of a circular plate held firmly at the edge 

[14] . It was desired that each psi of differential pressure produce at 
least 10 microinches per inch of strain. Transforming the equation of 
stress at the center of the — plate to the equation of strain (end 
of this appendix) and substituting the appropriate values, the plate thick- 
ness was calculated as being 0.0336 inches times the radius of the plate. 
One A-7 strain gage was located at the center of each side of the plate 
and the cumulative output of the two gages was read on a Baldwin strain 
indicator. This cumulative output was thus expected to be approximately 
20 microinches per inch of strain per psi of differential pressure. 

The plate thickness actually used was 0.032 inches due to the 
availability of this uniform thickness sheet stock. This plate was 
welded between two ertoW cylinderical shells of 1.5 inches inside 
diameter and 2.25 inches outside diameter. While this combination of 
plate thickness and effective plate diameter exceeded the calculated 
plate thickness to plate diameter ratio, the combined output still ex- 
ceeded the desired minimum of 10 microinches. The open ends of the 
cylinderical shells were threaded to receive a plug which contained a 


pressure tap. An O-ring seal was used between the plug and the shell 
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of the transducer. 

One leg of each strain gage was led through the shell wall and 
the drilled hole secured with epoxy cement. The other leg of the gage 
was soldered to the inside of the shell which served as a common ground 
between the two gages. The two protruding wires from the gages and a 
common ground wire, soldered to the outside of the shell, were then 
connected to the Baldwin strain gage indicator. 

The completed transducer was calibrated on a hydraulic dead weight 
tester. One end of the transducer was vented to the atmosphere and the 
other cavity filled with oil and connected to the dead weight tester. 

A pressure differential eerie from 0 to 100 psi, in increments of 
10 pounds, was then imposed with strains recorded while loading and un- 
loading the tester. 

A linear plot of pressure differential vs. strain was obtained for 
several transducers but on others, the welding caused a slight permanent 
buckle in the plate and the strain indications were not reproduceable 
for a given pressure differential. Due to the time involved in making 
and testing each transducer, it was felt that too much unproductive time 
would be spent on making transducers that could not be used; so this 
means of pressure ade butonent was abandoned in favor of manometer 
pressure measurements. 

One test transducer was glued together with Eastman 910 glue to 
eliminate the problems caused by the buckling of the diaphragm due to 


welding. A centering hole and shoulder was made in the cylinderical 


shells, the diaphragm was then glued to both shells under light pressure. 
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The strain gages were attached to each side of the diaphragm prior to 
joining the two end shells. This transducer gave a linear plot of 18.52 
microinches per inch of strain/psi of pressure differential when tested 
on the hydraulic dead weight tester but it did not respond to less than 
3/4 of an inch of mercury pressure differential when checked against a 
manometer reading of the same pressure differential with the tunnel in 
operation. 


The stress distribution in a uniformly loaded diaphragm with clamp- 


ed edges is given in Figure 1-D below [14] : 


DIAPHRAGM 
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Figure 1-D Stress Distributions of Diaphragm 
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The tangential stress at the center of the plate is given by the 
equation 

.. Seer) Pes 7 

eo 

Where S = Stress in psi. 
ie = Poisson ratio. 
'* = Pressure differential in psi. 
{’ = Radius of diaphragm in inches. 


t 


The corresponding value of strain is obtained from 


' Diaphragm thickness in inches. 


SD > 
.- - /z (2) 


where € is the strain and E = modulus of rigidity of diaphragm 


material. Substituting equation (2) into equation (1) 


Utilizing the required minimum value of € = ior inch / p54 
apcevential ge = 0.3 and E = 30x 10° psi for the values of 


steel the ratio of the diaphragm thickness to its radius may be obtained. 


—= 


jo = 3l«r'(/-0. 3" } 
B Re Biteed O” wrongly S 


eee | 
£ = 0.0336 a 


For a diaphragm radius of 0.75 inches a plate thickness of 0.0252 


inches is required. 
57 
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Figure 2-D. Exterior View of Welded Pressure Transducer, 
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Figure 3-D. Interior View of Welded Pressure Transducer. 
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APPENDIX E 
GENERAL THEORY OF THE SCHLIEREN METHOD OF OPTICAL INVESTIGATION OF 
COMPRESSIBLE FLUID FLOW {2] . | 

Apart from the conventional methods of experimentally investigating 
flow patterns vy means of pressure and velocity surveys, compressible 
gas flows lend themselves particularly well to optical methods of in- 
vestigation. 

Fundamentally, the optical methods in common use (the interfero- 
meter, the Schlieren, and the shadowgraph) depend on one of the two 
physical phenomena: (1) the speed of light depends on the index of 
refraction of the medium ‘PONE: which it passes, and the index of re- 
fraction of a gas in turn depends upon its density; and, as a conse- 
quence of this first phenomenon, (2) light passing through a density 
gradient in a gas (and therefore through a gradient in index of re- 
fraction) is deflected in the same manner as though it were passing 
through a prism. In a high-speed gas flow the density changes are 
sufficiently large to make these phenomena sizeable enough for optical 
observation. 

The interferometer, based on phenomenon (1), measures directly 
changes of density, and is primarily suited for quantitative determina- 
tion of the density field. 

The Schlieren method, based on phenomenon (2), measures density 
gradients. Although it is theoretically adaptable to quantitative use, 
it is inferior to the interferometer in this respect, and its greatest 


utility is in giving an easily interpretable picture of the flow field 
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together with a rough picture of the density variations in the flow, 

The shadowgraph method, also based on phenomenon (2), measures 
the second derivative of the density (i.e., the first derivative of 
the density gradient). Therefore it makes visible only those parts 
of the flow where the density, gradients change very rapidly, and it 
has found its greatest utility in the study of shock waves. 

Of the three methods mentioned, the interferometer yields the 
most information, aid the shadowgraph the least. On the other hand, 
the interferometer is the most costly and the most difficult to operate, 
whereas the shadowgraph is the least costly and the easiest to operate. 
Each method therefore has its own useful niche in experimental work, 
and the choice of method depends on the nature of the investigation. 


Figure 1-E below shows one of the numerous type of Schlieren 


(striae) arrangements. Light from a.uniformly illuminated line source 





TWO DIMEN SIONAL 
AOE L. = 






IMAGE OF 
a wg ay 












| FINITE 


SGREE N 
LINE SOURCE EDGE 


a: 


Figure 1-E Schematic Schlieren System 


60 


. ill 

a ll ey ea me 
ee ee ee 
ee ee 


= ee ae 6" be = « ee eh 


Ss _..  £ ies ae Ge i) Ga ae ae | 
hh, Settee a) Oe —eE——, 
oe — » mae aa, *' yy A 


———— © lS Ee Se CL es 5 az 





of small but finite width is collimated by the first lens and then 
passes through the test section. It is then brought to a focus by 

the second lens and projected on the screen, At the focal point, 

where there exists an image of the source, there is introduced a knife- 
edge which cuts off part of the light. With no flow in the test sec~ 
tion the knife-edge is usually adjusted so as to intercept about half 
the light, and the screen is uniformly illuminated by the portion of 
the light escaping the knife-edge. When the flow is established in 

the test section assumed, for simplicity, to be two imensional, with 
wach light ray pessing through a path of constant air density any light 
Lay passing through a regicn in which there is a density yradient nor- 
mal to the light direction will be deflected as though it has passed 
through a prism. Therefore, depending on the orientation of the knife- 
edge with respect to the density gradient, and on the sign of the den- 
sity gradient, more or less of the light passing through each part of 
the test section will escape the knife-edge and illuminate the screen. 

Thus the S:hlieren system makes density gradients visible in terms 
of intensity of illumination. A photographic plate at the viewing 
screen records density gradients in the test section as different shades 
of gray. 

For the Schlieren figure drawn it is assumed that the flow in the 
test section is parallel and in the x, y-plane, and that the light 
passes through the test section in the z-direction. 

Since the speed of a wave front of light varies inversely with 
the index of refraction n of the medium through which the light travels, 
it follows that a given wave front will rotate as it passes through a 


gradient in n. Accordingly, the normal to the wave front will follew 
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a curved path; since this normal is what we mean by the 
light ray, the latter is refracted as it passes through the 
density gradient. Noting that n is nearly unity for gases, 


it may be shown that 


2 gradient? 
BR 


where R is the radius of curvature of the light ray. The 
total angular deflection S of the ray in passing through 
the test section of width L is therefore given by 


> a = = l. grad /} 


Resolving this into Cartesian components, and connecting 
the index of refraction with the gas density through the 


empirical Gladstone-Dale equation which states that 


A! = | 
g Gap 
where Kop’ the Gladstone-Dale constant, is constant for a 


given gas, we get the angular deflections of the light in 


the x- and y-directions 


If the knife-edge is aligned in the y-direction (i.e., 
normal to the flow) as in Figure 1-E, only deflection in 
the x-direction will influence the light passing the knife- 


edge. Hence density gradients in the x-direction will be 
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made visible, but gradients in the y-direction will not be 
visible. Similarly, if the knife-edge is parallel to x, 
only gradients in the y-direction will become visible. 

In interpreting a Schlieren photograph of a two dimen-~- 
Sional flow, it is convenient to imagine that the photograph 
represents a view of a relief map of the flow field in which 
the vertical elevation is the gas density. If the map i's 
illuminated by glancing side light nearly parallel to x, the 
pattern of shadows and brightness will correspond to a 
Schlieren photograph with the knife-edge in the y-direction. 
Similarly, illumination in the y-direction corresponds to 
placing the Schlieren knife-edge parallel to the x-direction. 

Figures 14b and 14c are Schlieren photographs of the 
installed test section and show the two different effects 
obtainable with orientations.of the knitese@ge. In Figure 
14b (knife-edge horizontal) light areas on the upper half 
correspond to dark areas on the lower half, and vice versa; 
the boundary layers on the walls are clearly visible, since 
with this knife-edge orientation the density gradients nor= 
mal to the flow are made visible. In Figure 14c (knife edge 
vertical) the upper and lower halves look alike; the Houma 
ary layers are not visible, except aft of the shock wave 
which is inneracting with the boundary layer. 

The quality of focus of the Schlieren pictures is be- 
lieved to be poor due in part to the present arrangement of 
the Schlieren system installed for the reasons as outlined 
below by Holder and North [7]: 


Although it is usually desirable to arrange for the 
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beam of light passing through the working section to be 
parallel, nonparallel beams are occasionally used and an 
arrangement for producing this type is sketched in Figure 
2-E. Only one mirror or lens is required (apart from the _ 
focussing lens), but this advantage is to some extent off- 
set because its diameter must be larger in relation to the 


faetd to be examined than if a parallel beam is used. 
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Figure 2-E. Coincidence Schlieren System using 
a Plane Mirror or knife-edge 


In addition to employing a nonparallel light beam, 
the mirror system involves a double passage of the beam 
across the working section. The source is piaced at the 
center of curvature of the mirror, and part of the image 
of the source is reflected on to the viewing screen by a 


small plane mirror held near the source. Because the 
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source and its image are very close to each other, the 
arrangement is often referred to as the "coincidence" system. 
If the image of the source is deflected in a direction per- 
pendicular to the edge of the plane mirror, the illumination 
on the viewing screen changes; the edge of the mirror thus 
plays the same part as the knife edge in a Toepler appa- 
Batu’s . 

Since the source and its image are not truly coincident, 
the paths of the incident and reflected rays across the 
working section differ even in the absence of deflections 
produced by density gradients. This has an adverse effect 
on the image quality (in addition to any effect associated 
with the use of non-parallel light) since two slight offset 
images are produced. (See Figure 1l6c) 

These images may be made to coincide at the expense of 
considerably reduced illumination by using the arrangement 


sketched in Figure 3-E. 





Figure 3-E. Coincidence Schlieren System uSing 
Graded Filter 
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Here a graded filter is placed at 45 degrees to the axis at 
the center of curvature of the mirror. A condenser lens is 
used to produce an image of the source on the surface of 
the filter, and part of the light from this image is re- 
flected into the mirror. The image of the source formed by 
reflection from the mirror is made to coincide with the 
original image on the surface of the filter, and the light 
transmitted by the filter is arranged to fall on the view- 
ing screen in the usual way. The source and its image may 
thus be made exactly coincident so that, in the ahsence of 
density gradients, the zncident and reflected beams coin- 
cide. This will no longer be the case when the rays are 
deflected by density gradients, and some loss of image 
Sharpness, or difficulty in interpreting the observations, 
may then result. 

A further difficulty arises because only one of the 
two images formed by the incident and the reflected light 
can be in focus. It is usually best to focus the image 
associated with the reflected light since the incident rays 
are subject to further refraction during their second pass- 
age across the wind tunnel. 

The principal advantage of the system is that the 
double passage of the light beam, and the fact that the 
source image is formed at the center of curvature of the 
mirror instead of its focus, give increased sensitivity. 
For example, in a two-dimensional flow containing a 
region of density gradient of sufficient extent for the 


incident and reflected ray to remain within it, the 
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angular deflection of the ray will be doubled by the use of 
such a system. Moreover, the displacement of the image of 
the source for a given angular deflection wall be twice as 
large when the image is at the center of curvature as when 
it is at the focus. Coincidence systems may thus be val- 
uable where extreme sensitivity is required (for example, 
tunnels working with very low density), although for or- 
dinary use they are undesirable for the reasons mentioned 
above. 

Part of the increased sensitivity associated with a 
double-passage system may be retained, without the diffi- 
culties associated with the use of non parallel light by 
means of the arrangement sketched in Figure 4-E. Here a 
parallel beam of light is produced by a lens or mirror and 
after crossing the working section it is reflected back by 


means of “av plane mirror. 





tL OPLANE MIRROR 
SOURCE AND IMAGE AS IN FIGURE DE 


Figure 4-E. Schlieren System with Parallel Light Beams 
Crossing Tunnel Twice with Increased Sensi- 
Liviey, 
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angular deflection of the ray will be doubled by the use of 
such a system. Moreover, the displacement of the image of 
the source for a given angular deflection will be twice as 
large when the image is at the center of curvature as when 
it is at the focus. Coincidence systems may thus be val- 
uable where extreme sensitivity is required (for example, 
tunnels working with very low density), although for or- 
dinary use they are undesirable for the reasons mentioned 
aleve . 

Part of the increased sensitivity associated with a 
double-passage system may be retained, without the diffi- 
culties associated with the use of non parallel light by 
means of the arrangement sketched in Figure 4-E. Here a 
parallel beam of light is produced by a lens or mirror and 
after crossing the working section it is reflected back by 


means of a plane mirror. 





RR&A&Q@i\i OQ 


© PLANE MIRROR 
SOURCE AND IMAGE AS IN FIGURE Fé 


Figure 4-E. Schlieren System with Parallel Light Beams 
Crossing Tunnel Twice with Increased Sensi- 
C1V ey. 
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Tinsley Laboratories Inc. furnished the basic components for 
the coincidence Schlieren system as well as the design as 
shown in Figure 11. One unit was the six inch spherical 
mirror with a focal length of 80 inches mounted on a small 
pedestal. The other unit, mounted on a wooden frame with 
leveling jacks, consisted of the light source, knife edge, 
surface mirror and encased in a light tight wooden box the 
camera lens, folding surface mirrors and a six inch ground 
glass view plate. These were mounted on a portable stand 
as shown in Figure 9. 

An Ilex Optical Company Number 3 Acme Synchro Shutter 
was installed between the knife edge and the camera lens by 
means of a machined adapter fitted over the lens mounting 
with the shutter bolted to it. This arrangement allowed 
all of the light passing over the knife edge to enter the 
camera lens system but otherwise maintained the light tight 
integrity of the system. 

The six inch ground glass view plate was removed and 
the back of an Ilex Optical Company Korona camera consist- 
ing of a four inch by five inch ground glass view plate 
ama film pack holder fitted in its: places 

The procedure to align the Schlieren system is as 
follows referring to Figure 11: 

1. Locate the spherical mirror as close to the test 
section window as practical. 

2. Locate camera, light source, and knife-edge 


assembly such that the light source 1ie€S on an axis 
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perpendicular to the window. The light source should be 
approximately 80 inches from the spherical mirror surface. 
This may be done by moving the spherical mirror in the slots 
provided on the stand. 

3. With the light source "turned-On" adjust the 
spherical mirror to image the light source at the knife- 
edge. If the image is not in focus with the knife-edge 
assembly at approximately the middle of its slide travel. 
move the spherical mirror forward or back to correct the 
image focus. (A small black card is helpful in locating 
images and the light beam directions. ) 

4. Adjust the knife-edge screw such that the "Knife" 
does not intercept the light source image. The light should 
then pass through the camera lens, and an image of the test 
region should appear on the view screen. If the view screen 
image is not centered, adjust the small lst surface mirror 
and the spherical mirror such that the light passes through 
the center of the camera lens and aperture of the knife- 
edge assembly. If the view screen image is yet not centered, 
adjustment of the three small screw heads at the back of the 
camera will provide the final centering adjustment. 

5. The knife-edge can now be freely adjusted to provide 
the required Schlieren sensitivity. 

NOTE: All mirrors are front surface aluminized and 

should not be touched on the reflecting surface. 
They should be cleaned by dusting with a qual- 


ity camels hair brush. 
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The best photographs of the flow field were obtained 
uSing the following camera settings: 

1. Shutter speed set at 1/200 of a second. 

2. Lens opening at f-1l. This value has nothing to 

do with the actual exposure time of the film. 
This value of the lens opening allows all of the 
light beam to enter the camera and helps reduce 
stray light entering the lens system and fogging 
the film. 

3. Royal Pan film with a tungsten film speed rating 
of 160 was used. 

4. Developing time of eight minutes in total darkness 
inal: 1.dilute mixture of Eastman DK 50 devel- 
oper and water at 68°F. 

Colored Schlieren views of the flow field may be ob- 
tained by replacing the knife-edge assembly with the color 
grid purchased from Tinsley Laboratories. Colored pictures 
were not attempted as a means of printing and developing 
the color film was not available. 

The light source may be either a General Electric 
Company bulb number 1968 or number 1924, or a Westinghouse 
bulb number 1633. All three bulbs produce a line source 
of light. If the General Electric bulbs are used, a dif- 
ferent mount must be made for the bulb. 

Both the light source and the knife edge may be rotated 
from the horizontal position to the vertical position to 
obtain the desired direction of the density gradient of the 
flow field. 
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For the reasons already outlined, the focus of this 
Schlieren system is not as sharp as could be obtained with 
a two mirror one pass light beam system. 

A metal frame was fitted over one of the Schlieren 
windows and small copper wires were soldered into position 
at one inch intervals, starting at the throat position, to 


give horizontal dimensions. 
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APPENDIX F 
SAMPLE CALCULATIONS FOR EXPERIMENTAL DATA AND THEORETICAL 


PREDICTIONS. 


PREDICTION OF ADIABATIC MACH NUMBERS FOR NOZZLE AND CON- 
STANT AREA TEST SECTION OPERATING WITHOUT SHOCK. 

Flow is assumed to be isentropic to the nozzle throat 
and adiabatic thereafter. For ease of calculations, the 
nozzle contour as given in Figure 19a was assumed to con- 
Sist of five straight line segments as listed in Table X. 

Utilizing the influence coefficients of table 8.2 (2) 
multiplied by ue and designating the product F as outlined 


by Shapiro (2) for the combined effects of area change and 


friction results in the following equation 


w= 


2 DA 
= 4 
qdM = . ae FR4t 


O: 


ax. 
D 

The technique of isoclines as given in Chapter 8 of Shapiro 

(2) was the means of arriving at a solution to this differ- 

ome al seqwation. 

The parameters of the area A and the hydraulic dia- 
meter D may be expressed in terms of the tunnel dimensions. 
Area of the tunnel at any point is the product of the tunnel 
width w and the tunnel height h at the point in question. 

As the tunnel is of constant width, the area at any point is 
proportional to the height at that point and the derivative 
of the area proportional to the derivative of the tunnel 
height. The hydraulic diameter is defined as 4 times the 


area divided by the wetted perimeter and expressed in 
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tunnel dimensions is 


D = 4 x area - 4wh ge DN TE ae 
wetted perimeter 2 (w+ h) wth 


As the tunnel width w is constant throughout the length of 
the tunnel at 0.500 inches the hydraulic diameter may be 


expressed 
~~ 2(0.500) h _ 2 he. 
° t “ 1+2h 


As assumed, the nozzle segment heights will vary lin- 





early in accordance with the equation for a straight line 
h = mx+h, where mis the slope of the nozzle segment 
at some distance x downstream of the nozzle throat and the 


nozzle throat having an initial height of ho: 


h = MX + Ay 
dh= mdx > dx = be 
m 
hy, ~ Ay 
tox = LL, h = hy so the slope m = Tae 


Combining these relationships and substituting them into 


the original differential equation 


2 dh Ae L ( w+h dh 
= de goet Fe“! a ee 
F F_(4f) L 
h 2 Ww (hy 7 ho) h 
F.(4£)L 
am7 = *a 2 hai — i, (19s 
dh h h 


Z w (h,- h ) 


The equation is now in the form so that by suitable numer- 


ical values of h and M being substituted cm slope of the 
€ 
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Mach number squared is obtained. Values of h, are assumed 
in addition to values of M that are thought to bracket the 
actual Mach number. Ona plot of M2 vs. h, the calculated 
valves of the slope am*/dah are plotted at the coordinate 
points of each assumed value of Mand h. Sufficiently small 
increments of h must be utilized to provide the necessary 
isoclinic field. A smooth curve, starting tangent to be 
the valve of am*/dh at the known value of he and its corres- 
ponding value of Mis then drawn through the field of iso- 
clines tangent to a reference isocline at each value of h 
the curve passes through. In addition to the physical data 
listed in the previous table a value of the average friction 
fact on (4£) equaling 0.01 was assumed with an initial Mach 
number of 1.00 existing in the nozzle throat. 

A value of mM? is determined from the junction of the 
smooth curve and the height at the segment end point. This 
calculated Mach number is then used as the initial Mach 
number for the second nozzle segment. Continuation of this 
process resulted in the adiabatic Mach number at the noz- 
zle exit. Sample calculations for segment 2 are given in 
tabular form in Table XII and are plotted on figure 20. 

Straightforward use of Table 42(3) for Fanno line 
procedure was used to predict the resultant Mach number at 
the end of the constant area section of the tunnel. Iso=- 
clinic procedures could have been used but this method 
was discarded as being time consuming compared to the use 


of the tables. The following information was assumed or 
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utilized for determining the exit Mach number: o/ /[«.1 


aE = Q.01 and constant in value throughout the 
section 

I. = 5.75 inches, the tunnel section length 

D = Q.775 inches, the hydraulic diameter 

M) = wach-—number at the»enozzle~exit equal 2.017 

M. = Mach number at the section exit 


Entering Table 42(3) with M, of 2.17, a value of 4£L ay/P 


was read as 0.35298 for use in the following equation: 
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Me f f 2 
=> 0.35298 - (0.01) (5.75) _ 0.27878 
0.775 | 


Table 42(3) was then entered using this value of 4fL/D and 
the Mach number at the end of the test section determined 


as 1.914. 


PREDICTION OF ADIABATIC MACH NUMBER WITH A NORMAL SHOCK 


STANDING IN THE TUNNEL NOZZLE. 


Flow was assumed to be adiabatic in the nozzle until 
the normal shock occurred at a point midway between pres- 
Sure taps 4 and 5. M. the Mach number upstream of the 
normal shock was determined as 2.0 as read from figure 20, 
which is a plot of the adiabatic Mach numbers as calculated 
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vs. tunnel length based on previous calculations in this 


appendix. Entering Table 48(3) with this value of M the 


Mach number M, downstream of the normal shock was obtained 
Sore 77s5. M, thus obtained was used as the initial value 
of M for the isocline technique for determining Mach numbers 
for points downstream of the normal shock in the nozzle and 
in thewconstant area section of the tunnel. The nozzle area 
and the corresponding nozzle height were read from figure 
19a at the shock location assumed to be between pressure 
taps 4 and 5. Nozzle height at the assumed shock location 
is 1.30 inches. 7 

As before, the remaining portion downstream of the shock 
in the diverging nozzle was assumed to consist of the same 
srmpemght line segments. The initial section was taken as 
the remaining portion of the original segment 2 and was 
designated as segment 2a. All previous assumptions con- 
cerning the average friction factor and the physical data 
remained constant as outlined in Table XI. 

A nozzle exit Mach number was found to be 0.429. The 
relationship of 4fL/D at the inlet and exit of the constant 
area section was utilized to determine the Mach number at 
the end of the constant area section as was done in pre- 
vious calculations in this appendix. This exit Mach number 
was found to be M = 0.434. 


Theoretical values of P/P | at any point downstream of 


1 
the shock are obtained by multiplying the value of PAP 


corresponding to the Mach number at that point by the value 
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of eee obtained from Table 48 (3) for the Mach number 
at which the shock occurred. Based on the data in Table 
IX, it was assumed the shock occurred at a Mach number of 
2.000 with Bel og having a value of 0.7208. Adiabatic 
Mach numbers and pressure ratios obtained by this procedure 
are presented graphically in Figures 22 and 23. 

Table XVI is a Summary of the adiabatic Mach numbers 
predicted for the tunnel operating with and without a nor- 
mal shock standing in the tunnel nozzle with a constant 
area Seetion. 


TABLE XVI - Summary of Adiabatic Mach Numbers 
in Nozzle 


Pressure tap Adiabatic Mach Adiabatic Mach 
Number Number 
No Shock Shock in Nozzle 
2 Leo lao 
5 2a07 0 Om5 29 
7 2.150 0.432 
9 2. 195 0.429 
12 2.1m 0.429 
oi 1.904 0.434 


(end of constant 
area section) 
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DETERMINATION OF DIFFUSER HEIGHT TO PRODUCE CONSTANT MACH 
NUMBER IN TEST SECTION INCLUDING EFFECT OF FRICTION. | 

Assume that the diffuser height will increase linearly. 
Point t is taken as the nozzle throat; point 2, the end of 
the nozzle; and point 3 the maximum diffuser height loca- 
tion. . 


Utilizing the influence coefficients for constant 


specific heat and molecular weight from Shapiro's Table 


eee ie) 
ane 24 S41) dn, kevin FMI ELE 
Me =e A /=poae D 


DA. Zep? ai 
eo = AfkR/] D 


Area of tunnel = width x height = wh and as tunnel is of 


constant width 


Bin tt and % ~ dh 


/) 
ee ee eC a os Fee Pe 
ene te tamewer Wetted perimeter 2 (w+ h) 
2wh. Substituting the values 
w+ w+eh 


78 


__—e A a om 
— —— oe wat @ = Oars , 


ia 

— ee — a 
— +e =! = > —- ve 
in = 


2S 
~~ 





Combining terms 


4awdh _ gek idx 


W+h 
_— 2,7°*3 
FIT ( dx 


| 


os w+h 


a between designated limits 


h, - 
, (w+h) J, = BNF (X,- 4) 


[In ( #t4s ) = RIF (kaa) 
if OU 


clearing and solving for h, 


BM *f (Ka~Xz) 


h, =-W+ (Wth)e@ i 


2 


Substituting values of w = 0.5 inches; hoy = lovis ineh; 
M5 = Pawel s X3 - Xo = 5.33 inches, Figure 20. 


LE#2. re (0, 0025 NS,33} 


25 


0.5 + (2.218) © 


hz = 0.542.645 = 2.145 inches. 


The maximum theoretical height of the diffuser to ob- 
tain a constant adiabatic Mach number of 2.17 could not be 
obtained in the diffuser section due to physical limita- 


tions of the tunnel itself. As may be noted in Figure 19b, 
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the diffuser contours actually existing are far from the 
theoretical solution of the maximum. diffuser height. The 
actual contour required is of an exponential nature as 
developed previously. Values required in the diffuser to 
produce a constant adiabatic Mach number in the diffuser 
vs. diffuser length are presented in Table XIII in tabular 
form and plotted in Figure 19b. Numerical values of x were 
inserted in the exponential equation derived for determin- 


ing the diffuser height. 


DETERMINATION OF MINIMUM DIFFUSER AREA TO ALLOW SUPERSONIC 
FLOW. 

As in the case of the diverging nozzle, it was assumed 
the diffuser converged linearly from the maximum nozzle exit 
height of 1.718 inches to some minimum height less than 
1.718 inches. The same isoclinic technique was employed to 
determine this minimum height and allow supersonic flow to 
exist. The most severe condition that may be encountered 
is to have a normal shock at the nozzle exit or at some 
other location upstream of the diffuser throat. This con- 
dition will exist upon the initial start of the tunnel blow. 
The minimum diffuser height must be sufficiently large to 
pass the given mass rate of flow in the subsonic condition. 
Assuming a normal shock occurring at the nozzle exit, with 
a Mach number of 2.17 existing in the nozzle exit, the 
initial Mach number for the diffuser will be 0.551 as taken 
from the Gas Tables (3) for normal isentropic shock. At 


the minimum area, the Mach number must increase to 1.00. 
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Various values of the minimum diffuser height h are assumed, 
the isoclinic field drawn for these conditions and the 
smooth curve faired through the isoclinic field. Various 
values of h, the minimum diffuser height, had to be assumed 
to meet the end point conditions, namely, the curve had to 
have an infinite slope at the junction of h and the Mach 
number squared equaling one. Table XIV presents the data 
for determining the minimum height in tabular form and 
Figure 2~F is a plot of the isoclinic field for a minimum 
diffuser height of 1.40 inches. All isoclines are not 
plotted to reduce the clutter near the terminal point of the 
curves. The slope of dmM*/ah changes radically near the ter- 
minal point requiring a large number of isoclines to be 


plotted. 


DETERMINATION OF MACH NUMBER AT DIFFUSER THROAT WITH DIF- 
FUSER HEIGHT = 1.40 INCHES CONSIDERING FRICTION 

Conditions at nozzle throat are designated by t, at 
the nozzle exit by 2, and at the diffuser throat by 3. 

Utilizing the influence coefficients of Table 8.2 
multiplied by M* and designating the product F as outlined 
by Shapiro [2] for the combined effects of area change and 
friction 


2 tec, 
dM =h Go + ie at 


Area of tunnel = width x height = wh 
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“AS AX h_ for constant width tunnel 


dM? =f 2) + 6 af dt 


hydraulic diameter D = _4 x area =.= _4 wh = 
Wetted perimeter 2(w¢ h) 
2 _wh 
w+ h 





Substituting this value of I/D 


pa Lf ds. iC of. 
AM =f ol edh Or + 6. AF dr 


Assume that the diffuser height will vary linearly in 


accordance with the relationship 

h= mn th, 
where m is the slope of the diffuser and x is the distance 
down the tunnel. 


At 


X= SS = 4.802 inches, h = ho = 1.718 inches, 


and at 


x = 10.8 inches, h = h. - 1.40 inches 


m= hs — —he 7 L4-117/8 _ ~ 0 0530 


X,- Xz (0.8 -4.8 


dh = tndx “ak = ahh 


, = Ws pL Fe 
[B+ ae + SE] ah 
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and 


My = 2 ioe 


for determination of F 


i ; 
A = oa = /, 5549 inches 


f = 0 0025 


N71 - 2./7 


from Figure 20 


aa | 
Ma, =[f Fy + Se (ap + N(hy- hy) 


4.(D.0025) 
NM; ie rl ae In + - 0, 0530 sie OF (4-17) 


. 2 4 a 
M1. ere, = ~Zecre™ 1 Z2OTIC Fe 


Assuming a value of =. to determine M and the values of 


the influence coefficients from Table B-6 [2], a process 
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> iene ae = 


of iteration may be made to find the correct value of M, 


s LO3¢ a7 - {9 


Try fon 163 5 14 - a 
ea 7695 
cr - -—/2,038 
# 
Mo-M; = ~ O22 Fy p 0.01 92 F 


Mf = 470- 0 22/ (4.1635) — 0.0792 (12.038) 


mf = t:70 
I 2 sas /, 042 
= 1 Mx OC 
Try Pig= Owes ad /, 705 
, - 4 7649 
Et el PPOs 


4 
amy 


Ma = 4.10 - 0,221 (4.1649) - 0,07 92((2. 106) 


Mz, = 269 
Ms, 164 
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In addition to the analytical solution for determining 
the Mach number at the diffuser throat, an isoclinic field 
plot was drawn and an exit Mach number obtained by this 
technique. The Mach number as determined by the use of 
isoclines was 1.71 compared to a Mach number of 1.64 as 
determined by the analytical solution. Table XV and Figure 
3-F present the information necessary to plot the isoclinic 


field and a plot of the field itself. 


DETERMINATION OF EXPERIMENTAL ADIABATIC MACH NUMBERS FOR 
SHOCK FREE CONDITIONS ; 

Assume that flow to the nozzle throat is isentropic and 
that the flow for the remainder of the tunnel system is 
adiabatic. From the Gas Tables for isentropic flow (3) the 


adiabatic Mach number, M, at any desired point may be deter- 


mined for shock free flow conditions under the product 


(A) 2) | 


providing A* and P, are taken as the area and stagnation 


pressure when the Mach number is unity. P/P. = 
S -@e 
P 


O 
use of manometers. Consider pressure tap 16 at the outlet 


with pressures Py and AP being measured by the 


of the nozzle as station 1 and tap 31 in the diffuser as 
swativer 2. 

Data is taken from Table V or Figure 20 for a shock 
free run with the diffuser set at a constant height of 


1.718 inches. 
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(A/AL) om 2084 
(P/Po), = 0.0844 
(A/AZ) (P/Pg), = 2.34 x .0844 = 0.1975 
From the Gas Tables for isentropic flow (3) using the above 


value for entering the tables M, is determined as being 


1 
2.125. All adiabatic Mach numbers upstream of a shock were 


determined in this manner. 


DETERMINATION OF EXPERIMENTAL ADIABATIC MACH NUMBERS FOR 
SHOCK IN NOZZLE OPERATION. 

The adiabatic Mach number for points downstream of the 
shock in the nozzle sak ae eemeed in the same manner as 
for shock free conditions. Taking data from Table. IX for 
pressure tap number 20, the product (AJA) (B/Bo), was cal- 
culated as being 0.941. Entering Table 30 (3) with this 
value, the adiabatic Mach number was determined as being 
0.594. The ratio of change in stagnation pressure at pres- 


sure tap 20 or at any other location downstream of the shock 


in the tunnel to the initial stagnation pressure is obtained 


Fa) (9, 


being evaluated from table 30 (3) at 


by the relation 


=) 
Pol n 


The value of (P/Poy) 1 


the Mach number M, under the column P/P,. 
{ 
DETERMINATION ' OF EXPERIMENTAL MEAN FRICTION COEFFICIENT 
To eliminate as many variables as possible, data with 
the diffuser set at a constant height of 1.718 inches and 


the tunnel running shock free was used, i.e., Fanno Flow. 
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With the value of M having been determined from data in Table 
Vas being 2.125, at point 1, the constant area inlet, this 
Mach number was used as the argument to enter the Gas Tables 


(3) for Fanno Flow and evaluate the following: 


(=| = Ono fom 
pms) — Oa eae 
DD | 





BR _ (PIP), 
P (P/P") 


j 


we obtain 


- _ Of088 7, 2: 
am = SSS (0.3737) 
C 





0.48! 


\\ 


From the Gas Tables (3) for the Fanno Line flow using 
cay, - 0.481 the Mach number at the constant area 


pressure tap 31, point 2, is determined as being 1.78 and 
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— 


the corresponding value oe| pat as 0.23520 
Zz 


D 
el at) -0.3731- 0.2350 = 0.1384 
/ D Z 
m= O.7 75 inches 
Ai = 10.00. ~ 4.75 +=4.25 inches 
ae 
- 0 0/865 
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APPENDIX G 


AIR RECEIVER CAPACITANCE DESIGN CALCULATIONS AND PREDICTED 
TEMPERATURE VARIATIONS 


To maintain an isothermal blowdown, Muto [15] states 





Ta 
the ratio (I), ©.” = aoe should have as a minimum 
Pe Gy 
value of 7.5, and the non dimensional number of transfer 
units (2), NTU = LA. should have a minimum value of 7.5. 
‘Y) 
Vv 


Storage tank dimensions 


48 in, Outside Diam 
iS it ae in peaghe 
.75 in, heads 
.781 in, shells 
Mass in tanks = NM. /bm. 


BP pen (Tot! + 2 rat ] 


14222) ; 
ae Tl *® AT2Z* O.T8/* 1H $ Tee 60.75 | 
"STEEL (lee J 


M. = 10,410 16 mass of steel, 
Surface Area of tanks 
Z 
=2| 7d ‘ / ae 2A inde | = 
inside Z 


77 (4-6, 44)" 
- gp Tete ss Lge. a7) | 


/4.4 


AREA = /83.5 Square feef, 


25 





Mass of air initially in tank as beginning of blowdown = Mb 


Initial tank pressure = 475 psia. 


475'Vin®k lag Ye  fonks «177 Ypone 


Mb = = 
a 43.34 ft-# ,~% 620°R 
lhm -°R 
Mb = 577 Wem. adr initially in tank 


Mass of air in tanks at end of blowdown at 150 psia 


150 
= ——— = S/ 
Me = S77 —= = 182 sb, @ 150 PSIa 


MI a 
mM = 7 = "82 /577° tes 
ap 
Substituting the above values into equations 1 and 2 


- 








a a Me Cop sree _ 10,410 K OMF | 18 
ol . eal 
oor 577 * O'73 
uSing a value of h = 3 em and mM v [/00 lege 
NTU - A A = 2 (83.9 
Cy 0.173 xl/O00 


NTU = 2.89 
Since the value of ae 7? (7-5 and the value of NTU ¢ ies 
the value of NTU must be increased by increasing the area 


available for the heat transfer 


h A 
Cy m 


A = eae oui ‘* Diet 


NTU = 7.5 = 


94 





additional area required 
476 ff" ~ /83,5 ft" = 2925 ft 


if 1/8 inch aluminum sheet is used to reduce to corrosion 
problem, the mass of aluminum required utilizing both sides 


of sheet for heat transfer 


lim 


Q | Seats 2125 $4] = 169 FB [L525 $17} 





Mass Al 


Mass Al = 257 l1bm. 


Utilizing eqn 9A as developed by Muto [15] to predict the 
storage tank air temperature after a blow down with the 
presently installed tanks without additional corrective 
heating surfaces and comparing this value to the tempera- 
ture actually noted after a blowdown very close correla- 


tion was noted. 


Egn 9A by Muto to predict temp after blow 


- F ( Ate Gre . (231+ U1. B+ # 
~ 7G VF Co [+8 
T* = 0.9785 


with To = 523°R from Figure 1-G 
O 
Tend of blow = 0.9785 x 523°R = 511 R 
From Figure 1-G, actual temperature of the storage tank 


air at the end of a 10 minute run = 510.5°R. 
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ee 
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SF HE na FETE HH HEE Ce ce it ALE REA ETT 
a ae Fee ee eee ee eee eee TS 
SEER Mie HH SHS TTT He HH AEE Pet HEE Hea Hn 
ee HASLER eT Tt TURP LTE a THERES ESTEE Hl 
i 2) EERE LE AR SETHE 
/ HTT SE HELE TELS Aa a at IAEA EE i 
tate He THERA AGUAS UT AST ESTEE EST HE Hn HE HEE TREE HAE 
HET To SE ee ee CASITA i 
Hi ea PERE RED Pare LUE TMA Te Raed ea] ET 
FAME PR Hat Ce sree TEES AEE THE 
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